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Abstract The molecular phylogeny of the Pinaceae represented by nine species of eight 
genera was investigated through PCR-RFLP analysis of a cpDNA fragment, which was 
approximately 2550 bp long and homologous with the DNA fragment from the upstream of 
треї, to part of accD in cpDNA of Pinus thunbergii. Identified by 18 restriction 
endonucleases were 86 recognition sites, of which 54 were polymorphic. The restriction site 
data were analyzed with PAUP (version 3.1.1) and Mega (version 1.01). As a result, the 
most parsimonious Wagner tree constructed by branch-and-bound search is similar to the 
neighbor-joining tree in topology. Cathaya, Pinus, Larix and Pseudotsuga form a 
monophyletic group, and Larix is closely related to Pseudotsuga. Cathaya is more closely 
related to Pinus than to Pseudotsuga and Lariz, but the close relationship between 
Cathaya and Pinus is not well supported by the bootstrap test. Abzes has a close relationship 
to Keteleeria. Moreover, the system proposed by Zheng and Fu(1978) that the family 
Pinaceae was divided into three subfamilies, i.e., Abietoideae, Laricoideae and Pinoideae, is 
artificial, because the former two are not monophyletic. 
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摘要 AH PER 方法 分 别 从 松 科 8 属 \9 种 植物 中 扩 增 出 一 长 约 2550bp 的 cpDNA 片段 ， 这 一 片段 包 
括 rbcL\trnR、 部 分 accD 及 基因 间 的 非 编码 区 (相对 于 黑 松 cpDNA 中 的 同 源 片 段 而 言 )。 运 用 18 种 限 
制 性 内 切 酶 对 这 一 cpDNA 片段 进行 酶 切 分 析 ， 共 获得 86 个 酶 切 位 点 , 其 中 54 个 为 变异 位 点 。 运 用 
PAUP(version 3.1.1) 和 Mega(version 1.01) 软件 对 数据 进行 分 析 , 结果 Wagner 简约 树 和 Neighbor- 
Joining 树 反映 出 的 松 科 系统 发 育 关系 基本 一 致 : 银 杉 属 、 松 属 、 黄 杉 属 和 落叶 松 属 形成 一 个 单 系 群 ,， 且 
银 杉 属 与 松 属 的 亲缘 关系 更 近 于 与 另外 二 属 的 关系 ,但 这 一 结果 未 得 到 Bootstrap 分 析 的 较 强 支持 ; 落 
叶 松 属 近 缘 于 黄 杉 属 ;冷杉 属 近 缘 于 油 杉 属 。 此 外 , 松 科 中 的 冷杉 亚 科 和 落叶 松 亚 科 均 不 是 单 系 类 群 ， 
将 松 科 划分 为 冷杉 亚 科 、 落 叶 松 亚 科 和 松 亚 科 三 个 亚 科 的 系统 ( 郑 万 钧 , 傅 立 国 , 1978) 是 不 自然 的 。 
关键 词 ”分 子 系统 发 育 ; 限制 性 片段 长 度 多 态 性 ; rbcL 基因 ; 松 科 ; 银 杉 属 
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松 科 Pinaceae Lind. 是 现存 裸子 植物 中 最 大 的 一 个 类 群 ， 也 是 研究 松柏 类 植物 
(Conifers) 系统 发 育 的 一 个 关键 类 群 。 随 着 系统 学 的 深入 发 展 及 分 支 分 析 方法 的 引入， 
在 松 科 系 统 发 育 方面 已 积累 了 丰富 的 资料 : Harc (1987) 运 用 123 个 性 状 对 松柏 类 的 63 
个 属 (包括 松 科 所 有 10 个 属 ) 进 行 了 分 支 分 析 , 认为 松 科 是 一 个 很 好 的 单 系 类 和 群 , 种 子 中 
具 树脂 道 和 裂 生 多 胚 现 象 的 冷杉 属 Abies 、 金 钱 松 属 Pseudolarix ҢАР Keteleeria 、 雪 松 
B Cedrus 和 铁 杉 属 Tsuga 分 为 一 支 , 而 内 皮层 具 凯 氏 带 加 厚 且 次 生 木 质 部 中 其 树脂 道 的 
银 杉 属 Cathaya 、 松 属 Pinus 、 落 时 松 属 Lariz 、 黄 杉 属 Pseudotsuga 和 云 杉 属 Picea 形成 了 
另 一 支 。Price 等 (1987) 运 用 免疫 学 方法 对 松 科 9 属 植物 (不 含 银 杉 属 ) 进 行 了 系统 学 研 
究 , 并 将 松 科 分 为 两 个 群 :Abietoid group( 含 油 杉 属 \ 冷 杉 属 、 雪 松 属 、 铁 杉 属 和 金钱 松 属 ) 
和 Pinoid group( 含 云 杉 属 、 松 属 、 落 叶 松 属 和 黄 杉 属 )。 类 似 的 工作 虽 对 松 科 的 系统 发 育 
研究 有 一 定 的 启示 , 但 其 依据 的 性 状 存在 很 大 的 争论 。 松 科 在 地 质 历史 时 期 曾 发 生 十 分 
强烈 的 分 化 , 但 绝 大 部 分 类 群 现 已 绝 灭 , 现存 属 间 存 在 很 大 的 间断 , 具有 裸子 植物 中 普遍 
存在 的 孤立 性 , 很 难得 到 性 状 的 演化 系列 。 古 老 性 和 孤立 性 势必 导致 较 多 的 平行 演化 事 
件 及 性 状 的 相关 性 差 , 利用 形态 性 状 重建 其 系统 发 育 过 程 十 分 困难 。 八 十 年 代 以 来 , 分 子 
生物 学 技术 的 飞速 发 展 及 其 向 植物 分 类 学 这 门 古老 学 科 的 渗透 极 大 地 推动 了 植物 分 类 学 
的 发 展 , 利用 遗传 物质 (DNA) 本身 研 究 植物 间 的 系统 进化 关系 已 成 为 学 科 的 生长 点 之 
一 。 由 于 作为 DNA 基本 单位 的 碱 基 具 有 相对 稳定 的 变异 速率 ， 它 对 揭示 裸子 植物 等 古 
老 类 群 的 进化 过 程 起 着 尤为 重要 的 作用 。Chase 等 (1993) 利 用 rbcL. 基因 序列 对 499 种 种 
子 植物 进行 了 系统 发 育 重建 , 其 中 包括 松 科 9 属 (不 含 银 杉 属 )13 种 ;Tsumura $ (1995) 
对 松柏 类 ( 含 松 科 9 属 , 不 含 银 杉 属 ) 的 6 个 cpDNA 基因 进行 了 PCR-RFLP 分 析 ， 并 构建 
了 系统 树 。Chase 等 和 Tsumura 等 的 研究 均 基 于 高 度 保守 的 基因 区 ， 且 他 们 的 研究 中 均 
不 包括 特产 中 国 的 活化 石材 料 一 一 银 杉 属 。 本 文 运用 PCR-RFLP 方法 对 松 科 植物 的 
rbcL-accD 基因 片段 进行 分 析 , 探讨 松 科 植 物 的 属 间 关系 ,特别 是 银 杉 属 的 系统 位 置 。 


1 材料 和 方法 














表 1 材料 来 源 
Table 1 The origin of materials 
类 群 (缩写 ) 来 源 

Taxon( Abbreviation) Origin 
Abies holophylla Maxim. ( ABH) Botanic Garden, Institute of Botany( CAS) 
Pseudotsuga menziesii ( Mirbel) Franco (PSM) Botanic Garden, Institute of Botany( CAS) 
Picea smithiana (Wall. ) Boiss. (PIS) Botanic Garden, Institute of Botany( CAS) 
Keteleeria euelyniana Mast. (КЕЕ) Botanic Garden, Kunming Institute of Botany( CAS) 
Cathaya argyrophylla Chun et Kuang (CAA) Jinfo Mountain, Nanchuan, Sichuan 
Pinus armandi Franch. (PIA) Botanic Garden, Institute of Botany( CAS) 
Pinus banksiana Lamb. (PIB) Botanic Garden, Institute of Botany( CAS) 
Larix olgensis Henry (LAO) Botanic Garden, Institute of Botany( CAS) 
Cedrus deodara (Roxb.) G. Don (CED) Botanic Garden, Institute of Botany( CAS) 
Pseudolarix amabilis (Nelson) Rehd. (PSA) Botanic Garden, Kunming Institute of Botany( CAS) 


Tsuga chinensis (Franch. ) Pritz. (TSC) Botanic Garden, Kunming Institute of Botany ( CAS) 
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1.1 材料 来 源 

用 于 DNA 提取 的 材料 为 新 鲜 或 硅胶 干燥 的 针 叶 , 详 见 表 1。 
1.2 总 DNA 提取 

总 DNA 提取 方法 参照 Rogers & Bendich(1988) 的 方法 , 并 稍 作 改 进 ， 即将 1V CTAB 
沉淀 液 沉淀 DNA 这 一 步骤 改 为 TV 异 丙 醇 沉淀 DNA。 
1.3 PCR 反应 

扩 增 反应 在 PE 公司 的 9600 型 PCR 仪 上 进行 , 反应 体积 S0pl ， 内 会 50 mmol/L 
Tris-HCl (pH8.3), 0.25 g/L BSA, 2 mmol/L МЕСІ,, 1.5 u Taq DNA 聚合 酶 ( 购 自 北京 
农业 大 学 ), dATP, dGTP, dCTP 和 dTTP 各 200 pmol/L, 30—50 ng 总 DNA, 引物 PLA 
P; 各 12.5 pmols, P, 的 序列 同 黑 松 Pinus thunbergii rbcl 基因 5° 818 1 —25 TRÆ, 
Вр АТС TCA CCA ААА АСА GAG ACT AAA G °` (Wakasugi etal., 1994); P, 与 黑 松 
叶绿体 基因 组 41904 — 41925 区 域 (aacD 基因 中 ) 同 源 , 其 序列 为 ? ACT АСА GAT CTC 
АТА CTA CCC C? 。 扩 增 程序 为 :70 C ,4 min 一 94 C,2 min;50 C ,20 s;72 C,150 s; 
两 个 循环 ->94 С,20 s; 50 1С, 20 s;72 C , 150 s; 36 个 循环 ->72  ,6 min, 
1.4 ”限制 性 内 切 酶 消化 

实验 选用 Ava П. ВатнН І .Bgl II, Bsp1286 I 、BstN Í ,EcoR V .НаеШ. Hha I, Hind 
Ill Hint I „Msp I „Pst I „Rsa I „Sac I ,Sph I „Sty „Тад I „Xba T 35 18 种 限制 性 内 切 酶 
分 别 进 行 单 酶 切 ， 为 判断 一 些 疑 难 位 点 的 确切 位 置 ， 还 采用 了 BstN I + EcoR V Hha I 
+ BamH Í „Msp I + Bsp1286 I „Msp I + Hindlll „Rsa I + Sac I „Sty I + EcoR V Ж Taq I 
+ Hindl[ 36 7 种 双 酶 切 。 酶 解 反 应 体积 为 20pl, HA 3—5 p 1 PCR 397, 2—3 单位 限制 
性 内 切 酶 和 1 х 内 切 酶 反应 buffer。 酶 解 反 应 时 间 为 2 hr. o 
1.5 琼脂 糖 凝 胶 电 脉 

电泳 在 1.5% —2.0% 的 琼脂 糖 凝 胶 中 进行 ,电泳 buffer 为 1x TBE。 电 泳 结 束 后 , 直 
接 在 FOTODYNE 型 透射 紫外 灯 上 照相 。 
1.6 数据 分 析 

根据 限制 性 内 切 酶 位 点 的 有 无 得 到 1( 位 点 存在 )、0( 位 点 不 存在 ) 数 据 和 矩阵 , 然后 用 
PAUP 软件 (Swofford，1993 ) 中 的 分 支 界限 法 (branch-and-bound) 和 MEGA (Kumar et 
al., 1993) 软 件 中 的 Neighbor-Joining (NJ) 方 法 分 别 进行 分 析 。 运 用 PAUP 分 析 时 , 分 别 
用 变异 位 点 和 有 信息 的 位 点 (位 点 有 和 无 中 的 每 一 状态 至 少 被 两 个 以 上 分 类 群 所 共享 ) 参 
与 分 析 。 此 外 , 由 于 裸子 植物 科 间 的 间断 大 , 外 类 群 极 难 选择 (作者 曾 尝试 用 杉 科 植 物 作 
外 类 群 ,但 均 因 遗传 变异 太 大 , 参与 分 析 不 合适 ), 因而 PAUP 分 析 中 未 用 外 类 群 。 为 进 
一 步 判 断 Wagner 简约 树 中 各 分 支 的 可 信和 度 , 采用 了 Bootstrap 分 析 方 法 。MEGA 软件 分 
析 中 , 先 将 数据 的 1.0 状态 换 为 两 种 碱 基 状 态 , 再 以 P-distance 计算 遗传 距离 。 因 P- 
distance 计算 的 是 碱 基 替 代 的 百分率 , 因此 所 有 的 限制 性 内 切 酶 位 点 均 参与 分 析 。 


2 ”实验 结果 


由 Pi 和 Pz 扩 出 的 rbcl-accD 基因 片段 在 松 科 8 属 植物 (Tswuga 和 Pseudolarix 除外 ) 
中 无 明显 长 度 差异 , 均 长 约 2550 bp, 包括 rbcL 基因 1425 bp(Hipkins et al. ,1990), 其 余 
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表 2 86 个 限制 性 内 切 酶 位 点 在 松 科 8 属 \9 种 植物 中 的 分 布 
Table 2 The distribution of 86 restriction sites in nine species of eight genera in the Pinaceae 
_— Cy 





编号 num 。 。 识别 序列 Ша 变异 位 点 信息 位 点 分 布 (Distribution) 
NO. REP RS? СР? IS? АВН PSM PIS КЕЕ САА РІА PIB LAO CED 
character 

1 HinfI G! ANTC 36 1 dd doo d.d) d bd 
2 EcoRV GAT! ATC 106 £: ый, фо dq. dO Eo ch, do 71 
3 MspI C! CGG 138 0 0 0 0 1 0 1 0 9 
4 BsNI Cc! wGG 138 3 2 b; Л ¿de @ К, 3.00. 1-7 
5 HinfI G! ANTC 179 E do». dp $3.0. 21 
6 MspI C! CGG 189 £ Cy SS WS М do Q. Ck d 
7 Rsal GT! AC 193 p»UX docu OF SG d 4.3 
8 Аха G! GWCC 198 io do CE) od оо d. С 
9 Aval G! GWCC 210 pooh d dí Up 5d. eod 
10 Tagl T! CGA 231 4 3 0 0 0 0 0 1 1 0 0 
И TaqlI T! CGA 261 LU UR бй Тыс U cu t^ 
12 Haell GG! CC 266 5 4 0 1 0 0 0 0 0 1 0 
13 Bspl2861 | GDGCH! C 267 6 5 1.0 40 3 1 1! 0 0 14 
14 BstNI CC! WGG 273 7 0.4 n od od) uw d 
15 На G! АМТС 283 8 6 o 1 0 0 1 1 0 1 1 
16 Taql T! CGA 324 i; (£ qu ЛЫ» k UO E И 
17 HinfI G! ANTC 378 9 1.1 1 101 0 1» )! 1) 1 
18 Hael GG! CC 386 С da (b Po DER dod 
19 Bsp1286] GDGCH! C 397 10 0 1 1 1 1 1 1 1 1 
20 BamHI G! GATCC 417 11 7 0 1 0 0 0 0 0 1 0 
21 HinfI G! ANTC 418 12 8 1 0 1 1 1 ? 1 0 1 
22 Ауа G! GWCC 449 E qo. Ж” C od d Фу E 
23 Тад T! CGA 471 13 о о 0 0 0 0 о 1 0 
24 Haell GG! CC 498 14 9 1 1 1 0 0 0 0 1 1 
25 Rsal GT! AC 517 j as u OC Uu sb Rf 
26 Haell GG! CC 545 15 0 0 1 0 0 0 0 0 0 
27 HinfI G! ANTC 588 16 ] 3.1 14 X d 1.101 0 
28 Tal | T! CGA 591 17 1 ааа OE U O K 9 
29 Syl C! CWWGG 600 18 y- ае ob OO Op. 9 bod 
30 Hhal GCG! C 639 y LC BV d^ o4 Жуй 96 
31 PsI CTGCA! G 733 í do odo od do е d 
32 Rsal GT! AC 736 que S Tu АНЕ cb b o4 
33  Haell GG! CC 864 19 0 эф А7 do. hod 
34 HhaI GCG! C 887 1 Go do od) rp Pod: d. 21 
35 SphI GCATG! С 896 j' ur do cd ас ca? SE А 
36 Taql T! CGA 903 20 1 0 0 0 0 0 0 0 0 
37 Sphl GCATG! C 928 21 10 0 1 0 0 1 0 0 1 0 
38 Rsal GT! AC 939 L Gbed $^ ocE Cun ТИЕ ТИС. 
39 MspI C! CGG 963 p.d, OU Udo do Wu ds o1 X 
40 MspI C! CGG 984 22 n i dx ї do op o£ d d 
41 Rsal GT! AC 988 i- 40 X X0 AU odo ав 
42 TaqI T! CGA 1072 r Ды + $a cd de 11d 
43 BsNI CC! WGG 1116 £ 1 d Lo do do b dod 
44 Belll A! САТСТ — 1176 23 їз ДЕ 4. cd 6 1 lr 
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== Cont. 
ae num иин ma 变异 位 点 信息 位 点 M (Distribution) 

NO. RE” RS? СР” IS? АВН PSM PIS KEE САА PIA РІВ ГАО СЕР 

character 

45 На G! ANTC 1190 Беа о d.d: 1 
46 Rsal GT! AC 1197 24 p «codo do on £o oun ud 
47 Bsp12861 | GDGCH! C 1227 25 12 0 0 0 0 1 0 0 1 0 
48 SyI C! CWWGG 1229 їз Ges dp 240 р u a: GI 
49 Hhal GCG! C 1243 26 0 0 0 0 0 0 0 1 0 
50 MspI C! CGG 1245 27 0 1 0 0 0 0 0 0 O 
51 BstNI СС! WGG 1245 k Gc db Дф, ib und. d 
52 Taql T! CGA 1261 28 0 0 1 0 0 0 0 0 0 
53 Xbal T! CTAGA 1270 29 13 0 1 10 0 1 1 1 1 1 
54 Hindll A! AGCTT 1275 Lo d d È r do T a 
55 Rsal GT! AC 1284 I 1 id d Y Y T 4 1 
56 НіпаШ A! AGCTT 1341 30 14 1 0 0 0 10 1 1 0 0 
57 HinfI G! ANTC 1354 T dcr dod d d 1 1 
58 Ама G! GWCC 1455 31 0 0 0 0 0 0 0 0 1 
59 HinfI G! ANTC 1456 32 15 0 1 0 0 1 1 1 1 O0 
60 MspI C! CGG 1459 33 0 0 0 0 0 0 0 O0 1 
61 НаеШ GG! CC 1462 34 16 0 0 1 10) 1 1. 1 1. 1! 
62 HinfI G! ANTC 1475 35 17 1 0 1 0 00 0 0 0 
63  Haell GG! CC 1494 36 0 1 0 0 0 0 0 0 0 
64 ВЕШ A! GATCT 1510 37 18 0 1 10 1 10 0 0 1 1 
65 EcoRV GAT! ATC 1527 38 19 1- d qS 9 3. 3. 1. 0-1 
66 SwyI C! CWWGG 1600 39 20 0 0 0 0 0 1 1 0 O 
67 TagI T! CGA 1600 40 1 0 0 0 0 0 0 0 0 
68 TagI T! CGA 1651 41 1.1 10 10 14 0 3 1. 1 
69 SacI GAGCT! C 1710 42 B do or aic ЧУ Gd odo 1 20 
70 Bspl2861 | GDGCH! C 1720 1 1 10 10 10 1) 1. ! 1 
71 TaI T! CGA 1730 i de SPO Sp C О p. ka d 
72 Syl C! CWWGG 1740 43 o 0 0 0 0 0 0 O0 1 
73 MspI C! CGG 1753 i^ 3 ud ОП cw d Y 4d 
74 Aval G! GWCC 1755 p -ro do. do do de do X 
75 НІ GCG! C 1811 44 0 0 0 0 1 0 0 0 0 
76 BsNI CC! WGG 1860 45 0 0 0 0 1 0 0 0 O 
77 MspI C! CGG 1870 46 21 0 1 1 0 0 1 1 1 1 
78 BamHI G! GATCC — 1926 47 0 0 0 0 1 0 0 0 0 
79 Taql T! CGA 1990 48 22 1 0 0 о 0 0 0 0 1 
80 MspI с! CGG 2070 49 0 1 0 0 0 0 0 0 0 
81 MspI Ct cGG 2095 50 0 0 0 0 1 0 0 0 O 
82 Bspl286I — GDGCH! C 2170 51 23 0 0 0 0 0 1 1 0 O 
83 Таа] T! CGA 2190 £ 3^5 Xo do dd od dod 
84 Rsal GT! AC 2300 52 T God c0. do 1 Же od 
85 Xbal T! CTAGA 2340 53 24 1 1 1 0 1 1» 1 1 0 
86 Taql T! CGA 2375 54 25 1. 0 0 1 1 1 1.) 1 1 


1)RE: Restriction enzyme; 2) RS; Recognition site; 3)CP; Cutting position; 4) VS: Variable site; 5)IS: Informative site. 
长 度 包括 trnR、accD 的 一 部 分 及 基因 间 的 非 编码 区 。 此 外 , 各 属 的 PCR 扩 增 产物 中 均 有 
两 个 小 分 子 片段 , 一 个 长 约 1600 bp, 另 一 个 长 约 950 bp, 提高 扩 增 的 退火 温度 至 60 С, х 
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两 个 片段 也 不 消失 。Chase 等 (1993) 在 rbcL 基因 的 研究 中 发 现 很 多 类 群 具 有 rbcL 假 基 
因 ( 如 近乎 所 有 的 苏铁 类 植物 ), 松 科 中 的 这 两 个 小 片段 是 否 为 假 基 因 , 值得 进一步 研究 。 
Tsuga Fl Pseudolarix 较 特 殊 , 其 DNA 模板 在 其 它 基 因 ( 如 rDNA 中 的 26s 和 cpDNA 中 
的 trnk 基因 ) 的 扩 增 中 均 获 得 很 好 的 效果 ， 但 多 次 调节 扩 增 程序 均 扩 不 出 rbcL-accD 基因 
片段 , 估计 引物 Р, 不 匹配 ， 因 此 这 两 个 属 未 置 于 后 面 的 分 析 中 。 
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图 1 松 科 rbel.-accD 基因 片段 的 PCR-RFLP 分 析 1—8. Sac I ЮЛ; M. 100 bp DNA Ladder; 9—11. Hind Ш KHI. 
Fig.1  PCR-RFLP analysis of the rbcL-accD gene fragment іп the Pinaceae 1—8. SacI digestion; M. 100 bp DNA 
Ladder; 9 ~ 11. Hind WW digestion. 1. Abies holophylla; 2.9. Pseudotsuga menziesii; 3. Larix olgensiss 4. Picea 
smithiana; $. Pinus armandi; 6,10. Pinus banksiana; 7. Cedrus deodara š 8,11. Cathaya argyrophylla . 





图 2 松 科 rbeL-aceD 基因 片段 的 PCR-RFLP 分 析 1—8. Hinf I МИЛ; M. 100bp DNA Ladder; 9— 11. Hind ll НЛ. 
Fig.2 PCR-RFLP analysis of the rbcL-accD gene fragment in the Pinaceae 1—8. Hinfl digestion: M.100 bp DNA 
Ladder; 9—11. Hind Ш digestion. 1,9. Abies holophylla: 2. Pseudotsuga menziesii; 3, 10. Larix olgensis 4. Picea 
smithiana; 5. Pinus armandi; 6. Pinus banksiana; 7.11. Cedrus deodara s 8. Cathaya argyrophylla . 
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大 部 分 单 酶 切 结果 易于 分 析 , 如 图 1 及 图 2 中 的 9、10、11 三 个 样品 ;Hinf 工 的 酶 切 产 
物 中 含 很 多 小 分 子 片段 (图 2), 琼脂 糖 凝 胶 的 分 辩 率 不 够 ,结果 仅 分 析 了 rbcL-accD 基因 
片断 中 1 一 1500 bp 区 域内 的 位 点 ; 另 有 一 些 酶 的 切 点 不 好 判断 时 , 采用 了 双 酶 切 (图 3)。 
18 种 内 切 酶 的 单 酶 切 及 部 分 组 合 双 酶 切 共 获得 86 个 位 点 ， 其 中 54 个 为 变异 位 点 , 变异 
位 点 中 含 25 informative 位 点 。 所 有 内 切 酶 位 点 在 松 科 8 T8 .9 种 植物 中 的 分 布 见 表 2 
(Table 2)。 

















M1 234 5678 
. 





图 3 松 科 rbcL-accD 基因 片段 的 双 酶 切 

Fig.3 Double digestion of the rbcL-accD gene fragment in the Pinaceae М. 100 bp DNA Ladder; 
1,2. Larix olgensiss 3. Pinus armandi; 4~ 6. Pinus banksiana; 7—8. Cathaya argyrophylla . 
1,4. MspI; 2,5. MspI + Hind Ё; 3,6. Msp I + Bsp1286 I ; 7. Hha I ; 8. Hha I + BamHI 


用 于 PAUP 分 析 的 数据 矩阵 见 表 3。 用 54 个 变异 位 点 和 25 个 informative 位 点 分 别 


表 3 用 于 PAUP 分 析 的 数据 矩阵 
Table 3 Data matrix used in the PAUP analysis 





类 群 (Taxon) 性 状 (Characters) 
Abies holophylla 101001001001010111010111000001000010010111000001000111 
Pseudotsuga menziesii 101010111110010111101110001010010001110011000100100110 
Picea smithiana 001001101101011111100011000110000110110011000100000110 
Keteleeria evelyniana 100001101101000111100111000000000100100011000000000001 
Cathaya argyrophylla 110001110101000111101001100011010100110011011010010111 
Pinus armandi 101101111101000111100111000011020100011001000100001111 
Pinus banksiana 110100101101000110100111000011010100011011000100001111 
Lariz olgensis 101010111110110111101111110010010100100011000100000111 
Cedrus deodara 101001111101010001100111000010101100110010100101000101 


进行 分 析 , 各 得 到 一 个 最 简约 的 Wagner 树 ， 且 两 个 树 的 结构 完全 一 致 (图 4)。 运 用 
Bootstrap 检验 方法 对 数据 进行 100 次 重 抽样 分 析 , 结果 强烈 支持 松 属 二 种 间 的 亲缘 关系 
以 及 落叶 松 属 和 黄 杉 属 间 的 亲缘 关系 。 用 NJ 方法 构建 的 系统 树 如 图 5 所 示 , 与 Wagner 
简约 树 基本 一 致 


104 Bow 分 类 学 HW 35 卷 





93 Pseudotsuga menziesii 
Larix olgensis 
Cathaya argyrophylla 


Pinus armandi 





84 
Pinus banksiana 


Cedrus deodara 





Picea smithiana 








| Keteleeria evelyniana 





Abies holophylla 





图 4 ”基于 rbeL-accD 基因 片断 的 PCR-RFLP 分 析 , 运用 PAUP (version 3.1.1, Swofford, 1993) 的 分 支 界限 法 构建 的 
最 简约 Wagner Bl, 其 步 长 为 72, C1=0.750, H1=0.250, RI= 0.526, RC= 0.395。 分 枝 上 的 数值 示 利用 自助 法 分 析 


(Bootstrap analysis) 检验 某 一 分 枝 在 100 次 重 抽样 分 析 结 果 中 出 现 的 次 数 ( 仅 列 出 大 于 50 的 数值 )。 

Fig. 4 The most parsimonious Wagner tree generated by the branch-and-bound option of PAUP (version 3. 1. 1, 
Swofford, 1993) using PCR-RFLP analysis of the rbcL-accD gene fragment. The tree has 72 steps, with CI = 0.750, HI- 
0.250, RI= 0.526 and RC = 0. 395. Numbers above the branches indicate the percentage of times that the branch was 
recovered in 100 bootstrap replications. Bootstrap values below 50 are not reported. 
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= : this length is approximately equal to the genetic distance of 0.00272 


图 5 基于 rbeL-accD 基因 片断 的 PCR-RFLP 分 析 , 由 Neighbor-Joining 方法 构建 的 松 科 8 属 .9 种 植物 的 树 系 图 
Fig.5 The neighbor-joining tree of nine species of eight genera in the Pinaceae constructed from PCR-RFLP analysis of the 
rbcL-accD gene fragment. 
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3 讨论 
3.1 HR IR Sr dE 

本 研究 虽 仅 涉及 一 长 约 2550 bp 的 cpDNA 片断 ,但 86 AIRHEAD 8 54 
个 变异 位 点 ,已 占 63% 。 一 般 情 况 下 , 如 此 高 的 变异 量 对 揭示 系统 发 育 关系 极为 有 效 , 但 
本 文 分 析 的 54 个 变异 位 点 中 , 29 个 位 点 的 变异 为 自 征 , 仅 25 个 位 点 ( 占 变 异 位 点 的 
4696 ) 对 揭示 系统 发 育 关系 有 意义 , 这 无 疑 反 映 了 松 科 各 属 的 孤立 性 , 同时 也 是 Bootstrap 
值 偏 低 的 原因 之 一 。 
3.2 银 杉 属 的 系统 位 置 

银 杉 属 的 系统 位 置 一 直 存 在 很 大 的 争论 :其 胚胎 发 育 、 花 粉 形态 、 染 色 体 数目 和 核 型 
均 接近 于 松 属 , 木材 解剖 结果 接近 于 黄 杉 属 ( 王 伏 雄 ，1990)。 何 关 福 等 (1981) 通 过 化 学 
成 分 的 研究 , 认为 银 杉 属 近 缘 于 松 属 和 云 杉 属 , 这 一 结论 与 Hart(1987) 的 分 支 分 析 结果 
一 致 ;Chen 等 (1995) 对 银 杉 的 胚珠 结构 和 肉 配 子 体 发 育 进 行 了 研究 , 认为 银 杉 属 和 松 属 
具有 较 近 的 亲缘 关系 。 本 文 的 Wagner 简约 树 和 NJ 树 均 表明 : 黄 杉 属 落叶 松 属 、 银 杉 属 
和 松 属 构成 一 个 单 系 群 , 银 杉 属 与 松 属 的 亲缘 关系 更 近 于 与 另外 二 属 的 亲缘 关系 , 但 这 一 
结果 未 得 到 Bootstrap 分 析 的 较 强 支持 。 
3.3 落叶 松 属 和 黄 杉 属 . 冷 杉 属 和 油 杉 属 的 关系 

落叶 松 属 和 黄 杉 属 具有 极 近 的 亲缘 关系 (图 4 和 图 5), 冷杉 属 和 油 杉 属 具有 较 近 的 
亲缘 关系 (图 5)。rbcL 基因 序列 分 析 (Chase et al.，1993) 和 松柏 类 的 6 个 cpDNA 基因 
的 PCR-RFLP 分 析 (Tsumura, 1995) 也 得 出 这 一 结论 。 此 外 , Hart (1987) 的 分 支 分 析 和 
Price 等 (1987) 的 免疫 学 分 析 同 样 支持 上 述 结论 。 
3.4” 松 科 中 的 亚 科 划分 问题 

郑 万 钧 和 健 立 国 (1978) 将 松 科 划分 为 三 个 亚 科 一 一 冷杉 亚 科 、 落 叶 松 亚 科 和 和 松 亚 科 。 
在 这 一 分 类 系统 中 , 落叶 松 属 、 金 钱 松 属 和 雪松 属 属于 落叶 松 亚 科 ; 冷杉 亚 科 含油 杉 属 、 
冷杉 属 、 黄 杉 属 、 铁 杉 属 、 银 杉 属 和 云 杉 属 ; 松 亚 科 仅 含 松 属 。 冷 杉 亚 科 和 落叶 松 亚 科 的 主 
要 划分 依据 为 枝 是 否 分 为 长 短 枝 。 但 是 , Hart (1987) 的 分 支 分 析 发 现 短 枝 这 一 性 状 在 松 
科 中 至 少 起 源 了 三 次 :一 次 在 银 杉 属 、 松 属 和 落叶 松 属 中 ;一 次 在 雪松 属 中 ; 另 一 次 在 金钱 
松 属 中 。 短 枝 这 一 多 次 起 源 的 性 状 作为 亚 科 的 划分 标准 显然 是 不 合适 的 。 从 本 文 的 
Wagner 简约 树 和 NJ 树 上 看 , 落叶 松 属 和 雪松 属 分 别 位 于 不 同 的 支 系 中 , 冷杉 亚 科 也 不 是 
一 个 单 系 群 , Chase 等 (1993) 和 Tsumura 等 (1995) 构 建 的 分 支 图 也 不 支持 冷杉 亚 科 和 落 
叶 松 亚 科 为 单 系 群 , 说 明 郑 万 钧 和 傅 立 国 (1978) 将 松 科 分 为 冷杉 亚 科 、 落 叶 松 亚 科 和 松 亚 
科 的 分 类 系统 是 不 自然 的 。 
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